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We analyze the sector of dimension-three vector meson operators in the "hard wall" model of 
holographic QCD, including the vector and axial currents, dual to gauge fields in the bulk, and the 
tensor operator ipa^if), dual to a two-form field satisfying a complex self-duality condition. The 
model includes the effect of chiral symmetry breaking on vector mesons, that involves a coupling 
between the dual gauge field and the two-form field. We compute the leading logarithmic terms in the 
operator product expansion of two-point functions and the leading non-perturbative contribution to 
the tensor-vector correlator. The result is consistent with the operator product expansion of QCD. 
f*"^ We also study the spectrum of vector mesons numerically. 
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I. INTRODUCTION 



Popular phenomenological models of QCD, such as the "hard wall model" [l], HJ as well as "the soft wall model" jH, 
rely on the assumption that one can map QCD to an effective theory in the bulk of the holographic fifth dimension. 
O ,■ This is a very strong assumption, which is not fully justified. Such an effective description implies a large hierarchy of 
scales between meson masses and flux tube tension, which is not present in QCD. This is the main reason why these 
are at best phenomenological models. Nevertheless, they are a useful resource as they provide quick and easy ways 
to estimate many quantitative properties, provided one is willing to live with errors which, in most instances where 
these models can be compared to data, turn out to be in the 10-30% range. These simple tools are nice to have for 
QCD. They are even more valuable when studying QCD-like theories in particle physics, most notably as a theory 
of technicolor, or more generally as a potential "hidden sector" which may leave an imprint on LHC data. As far as, 
for example, meson spectra in QCD are concerned, holographic models will never be competitive with lattice gauge 
theories. However, when exploring theories of technicolor or hidden sectors one does not know, a priori, what the 
correct Lagrangian is. So one needs to explore many different models, each of which would require years of extensive 
computer simulations on the lattice, but only days in a holographic model. 

While using a 5D effective theory for QCD is not necessarily justified, it should at least be done self-consistently. In 
the 5D effective theory, not only are higher derivative interactions and high dimension operators suppressed — one is 
also suppressing an infinite number of additional fields. In holography, we know that every boundary operator should 
correspond to a field in the bulk. Fields kept in the simplest holographic models correspond to boundary operators 
of UV dimension[22[ 3. One can say that bulk fields dual to operators of higher dimension, which are more massive 
and hence can be integrated out, are being neglected. In a top-down holographic theory, such as the AdSs x S 5 dual 
of large Af = 4 SYM with a large number of colors N c and at strong 't Hooft coupling A > a similar reduction 
to a small subset of boundary operators and hence bulk fields is entirely justified: there is only a finite number of 
ten-dimensional fields dual to BPS operators, which retain their free field dimension. All other operators acquire 
anomalous dimensions of order A 1 / 4 , so their dual fields in the bulk acquire masses of the same order in A and can 
safely be integrated out. In the phenomenological approach one includes the fields dual to dimension 3 operators, but 
hopes to be able to neglect operators of dimensions 4, 5, 6 and so on. 

There are, however, two additional dimension 3 operators, 'ipa^ip and vbo^^^ (cr' Il/ = i/2[7'', 7"] being the 
antisymmetrized product of two gamma matrices), whose dual field is not included in the simplest bulk model. The 
corresponding field, a complex bifundamental anti-symmetric rank- two tensor field B^ v , should be included for self- 
consistency of the model. Real and imaginary part of this field correspond to the tensor operators with and without 
75 insertion respectively. An immediate benefit resulting from the inclusion of this extra field is that one will obtain 
masses for isospin triplet vector mesons with J PC = 1 H , starting with the 61 at a mass of 1235 MeV. They clearly 
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should be part of the setup which, as it stands, can otherwise only incorporate 1 and 1 ++ vector mesons like the 
p and Oi< 

Including a new held in the Lagrangian comes with new interaction terms and coupling constants. The original 
work on the hard wall model in ref. |l| has proposed a rigorous procedure to hx those, which so far has been very 
successful phenomenologically: first of all, we assume that the bulk is described by an effective field theory so we 
only write down interaction terms of bulk field theory dimension 5 or less (not to be confused with the dimension 
of boundary operators, which after all maps to the mass of the bulk field). The corresponding coupling constants as 
well as the normalization of the kinetic terms and the masses are obtained by demanding that at large momentum 
correlation functions in the bulk agree with the corresponding field theory values at weak coupling. This is basically 
the statement that one trusts the effective field theory picture in the bulk all the way into the UV. While this is a 
very strong assumption, it at least is a hypothesis that can be tested as it allows one to make predictions for particle 
masses and decay constants based on very few inputs. In the original hard wall model this procedure was used to 
fix the mass of the vectors and the scalar in the bulk, the five-dimensional gauge coupling and, as was pointed out 
later in ref. Q, it was also needed to fix the map between the asymptotic form of the scalar field and the quark 
mass, as well as the quark condensate (tptp) in the field theory. Out of these, only the mass of the vectors could have 
been justified 23] without the assumption of a valid effective theory in the UV. But already for the determination 
of the gauge coupling one had to rely on the two-point functions of the currents, whose overall normalization is not 
protected. 

Preliminary advances in including the B^ v field in the bulk has recently appeared in the pioneering work Q, but 
only in a form which does not account for chiral symmetry breaking. In that work only the real part of B^ v is 
considered; this field propagates the new l" 1 tensor mesons, but also an additional copy of 1 vector mesons. More 
importantly, only quadratic terms were included in the action for the new £?-field. Its mass is fixed to 1 in AdS 
units by requiring that the dual operator has dimension 3; the normalization of the kinetic term is fixed, as in the 
original hard wall model, by requiring the large momentum limit of two-point functions to agree with asymptotically 
free QCD. However there is one more bulk operator of dimension 5 or less that needs to be included in the action 
to communicate the effects of chiral symmetry breaking to the field and ensure that there is only a single set 
of vector mesons and no double counting, that would happen since gauge fields and the B^ v fields have the same 
degenerate spectra in the absence of this term. This extra term has the form Tr(X* F^B + BFrX^ + h.c), where X is 
the bifundamental scalar responsible for chiral symmetry breaking and F L / R are the field strengths for the bulk gauge 
fields dual to the chiral symmetry currents. The coupling constant in front of this term may be fixed by demanding 
the correct OPE structure of the correlator in the UV. 

More progress in this direction was made in ref. 0] , where it was proposed that the action of the complex £? M „ field 
should be first order, in such a way that the four-dimensional components satisfy a complex self-duality condition. The 
reason behind this choice is that in four dimensions the tensor operators tpa^tp and ipa^j^rp are not independent, 
but given the definition of 75 = i7°7 1 7 2 7 3 , they are related by 

^"75^ = ~/ V a ^ a ^- (1) 

Hence a similar condition must be imposed on the two-form field. Following a similar procedure, we choose the action 
for the B^ u to be a Chern-Simons action with a mass term.[24j schematically 

S B = - J d 5 x[ i{B A dB^ - B f A dB) + m B \B\ 2 ] . (2) 

The duality condition follows from the equations of motion. We differ from ref. @ in several ways. Just like in 
ref. Q, the authors do not include the effects of the dimension 5 bulk operator that communicates chiral symmetry 
breaking to the tensor sector, even though they correctly point out that its effect should be included. Secondly, in 
order to fix the degeneracies resulting from this truncation to a free field theory, they let the mass of the bulk B^ v 
field take different values, so that the field is dual to operators of dimension A different than 3. While it is true, as 
we pointed out above, that using matching to free field theory is not a well justified procedure, it is at least a testable 
assumption and so far has met with surprisingly large phenomenological success. If one abandons this, one should not 
just treat the mass of the B^ v field as a new free parameter, but also the five-dimensional gauge coupling, the mass 
and normalization of the bulk scalar field as well as the normalization of the B^ kinetic term, all of which affect the 
correlation functions in the boundary theory. In this case the model loses virtually all predictive power. Given the 
surprising accuracy with which the hard wall model so far has predicted particle masses, we believe it is premature to 
abandon the procedure of matching to UV correlators at this stage. We fix our parameters to reproduce the boundary 
expansion of a field dual to an operator of dimension A = 3. 

In this paper we will explicitly carry out the calculation for the short distance behavior of the bulk correlation 
functions of B^. Comparing to the operator product expansion (OPE) of weakly coupled QCD we will indeed be 
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able to completely fix all new coupling constants in the bulk. In fact, the set of conditions we obtain for the couplings 
is overdetermined and the fact that we can find values that allow us to reproduce all QCD correlation functions 
to leading order is a nice consistency check. The upshot is that this improved model has no new undetermined 
parameters. For now this serves as a proof of principle that this matching can be done. One has an improved hard 
wall model with no new free input parameters but several new predictions (masses and decay constants for the tensor 
and axial tensor mesons). Whether the phenomcnological success of the model survives these additions will be a 
good test to what extent the underlying assumption of an effective description in five dimensions gives an accurate 
picture of real QCD. We analyze the meson spectrum, we demonstrate that the cubic coupling indeed removes all the 
unwanted degeneracies of masses that were present in the case of a free B^, but unfortunately the meson spectrum 
that we observe does not match with what has been measured for QCD. 

The organization of this paper is as follows: in section|TI]we will review the short distance structure of the correlation 
functions involving the dimension 3 vector operators in QCD, as this is what we want to reproduce. In section Hill 
we present the improved holographic model. We derive equations of motion and the renormalized action in sections 
IIVI and [V] respectively. In section [V] we calculate the short distance correlation functions in the theory with massive 
quarks and extract the bulk coupling constants from comparing to QCD. As some of the correlation functions have 
the leading short distance terms proportional to the mass, in the chiral limit several correlators are dominated by 
the subleading term involving the chiral condensate. As our bulk Lagrangian at this stage is entirely determined, 
reproducing these correlators is a non-trivial check of our construction. We demonstrate that this indeed works out 
in section fVTl In section fVIII we analyze the meson spectrum and the phenomenological issues that appear in the new 
model. In section [Villi we summarize our results. 



II. CORRELATION FUNCTIONS IN QCD 



In two-flavor QCD, the relevant two-point functions in the vector sector are: 

n^y ab ( q 2 ) = % J d i xe iqx (n\T{V'" a (x)v v ^(o)}\n), (3) 

Il£7' a V) = * / d 4 xe iqx {n\T{T v P a (x)V^ h {0)}\n), (4) 

II^ Q/3 < a V) =* / d i xe lqx {n\T{T^ a {x)T a ^ b {Q)}\9), (5) 



where V fJ,a (x) = %l)(x)^^T a ip( x ) an d T^ va (x) = ■ijj(x)a ,J,v T a ip(x) are the vector and tensor isospin triplet currents 
respectively, and \Q) is the non-perturbative vacuum. We choose a normalization for the isospin generators such that 
tr (r a r b ) = \5 ab . The two-point functions above have the following kinematic structure: 

Iir/V) = 6 ab (q^-q 2 V n^Vv(q 2 ), (6) 

U^' ab (q 2 ) = S ab U+ T (q 2 )F^ afi +5 ab U^ T (q 2 )F^ a/3 , (7) 

n^°V) = i«$°VV - r,^q u )U VT (q 2 ), (8) 
where, defining the projector q 2 P tlv = q 2 rj^ v - q^q v , 

pa p P _ J_ pafi . / q s 

projects onto positive parity. Its counterpart is 

F^ af) = Fp afi - fW^rfP - ?y M V a ), (10) 

so F- ~ F+ — 1 is the negative parity projector. Notice that the sign of F_ is chosen so it is actually minus the 
projector 

FW°P = -{qVqPr,^ + q»q a lfP - q v q a if p - q^q^^ a ). (11) 

We also have: 
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with 



F±^ a8 F ± a i p = ±2q 2 F ± ^ p . (13) 



In the large- N c limit the two-point functions above are saturated by single-particle exchange of an infinite number of 
stable mesons, in this approximation to real QCD we can write, up to subtractions, the two-point functions above as: 

n n P> n 

T 



1 b 

with the decay constants defined as: 

J> , 



TT+ I 2\ \ ^ fb,n tt i 2\ \ fp.nfp n 



(n\V^\p b n (p, A)) = M p ^8 ab f p , n ^{ Pl A), (15) 



(n\T"\p b n (p, A)) = i6 ab fT n \p^e v (p, A) - p^( Pl A)], (16) 



(fi|2£,|6* (p, A)) = i5 ab f b , n e^ afi p a eP(p, A). (17) 

As it is made explicit by the notation above, the current V^ 1 produces vector mesons (J PC = 1 ) like the p, while 
the tensor operator produces both vector mesons and their even-parity partners (J PC = l" 1 ), like the &i meson. 
For large Euclidean momentum Q 2 = — q 2 — > oo contributions to these correlators can be organized according to the 
operator product expansion (OPE), with a leading perturbative contribution plus an expansion on the several vacuum 
condensates, (tpil)), (a s G 2 ), etc., that capture the non-perturbative effects. This was or iginally done for three colors 
in refs. |1(ML3 |. Expressions for a general number of colors can also be found in refs. [13l4l5l | . To leading order we 
have: 

q ^ nM<3 *) = -^,„ g ^ +0 (|l), (18) 

JSJ.^) (31) , , 19 , 



Q 2 ^oo lbn^ p z Q z \Q 



hm U VT (Q- 2 ) = --^mlog^-^+0[^), (20) 



where m is the quark mass. This is the large momentum behaviour of the correlators that we will use to fix the free 
parameters of the five-dimensional action. 



III. IMPROVED MODEL OF HOLOGRAPHIC QCD 

The model we consider is an extension of the hard wall model of ref. [l[, but it can be generalized to other 
holographic QCD models like the soft wall model of ref. p|. We use a five-dimensional geometry to describe the 
dynamics of four-dimensional QCD with a large number of colors N c — > 00. The metric is that of AdS^ with a radius 
£, we choose a mostly minus signature and work with the coordinate system 

£ 2 

ds 2 = g MN dx M dx N = — (~dz 2 + i^ v dx^dx v ) . (21) 

In these coordinates the boundary is at z — 0. In order to recover some of the physics of confinement, we introduce 
a cutoff in the radial coordinate z m . Since the radial coordinate maps to a renormalization group scale in the dual 
theory, with z — corresponding to the UV, the cutoff l/z m can be interpreted as an IR scale where the theory 
becomes confining. 

We introduce a set of fields cf>(x, z) in the five-dimensional theory that are dual to mesonic operators 0{x) in the 
field theory with conformal dimensions A < 3 and spin J < 1. In previous works this was done considering scalar 
and vector fields. This included both scalar and pseudoscalar mesons, as well as vector mesons 1 and axial vector 
mesons 1 ++ , but the full set of vector mesons include also l" 1 states, that were missing in the original formulation. 
These can be included by considering a complex two- index antisymmetric field B pu , or two- form for short. 
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The five-dimensional theory has a U(2)l x U(2)r gauge symmetry, that maps to the global flavor symmetry of 
two-flavor QCD. The fields A^^ and Ar^ will be the associated gauge bosons, while the complex fields X and 
are in a bifundamental representation 

X — ► U L XU R , — ► U L B^UI (22) 

The map between operators and fields can be summarized as: 

AD:Q(x) hD:<j){x,z) A m 2 / 2 

^l>tVi A a Lfl 3 

y2? fl 7,T°fe 3 

f L < 3 -3 

fo^i B$ 3 1 

Where m 2 JL 2 is the mass of the field. We have chosen masses such that the conformal dimension A of the dual operator 
matches with its free value. Although quantum corrections will change the conformal dimension of operators in the 
IR, QCD is a free theory in the UV and is in this regime where we will do the matching with our model, hence our 
choice of masses both for the scalar and the two-form field. We can also form the real combinations 

AD:0{x) 5D:<j)(x,z) A m 2 J 2 



^? = (A^+A£ /t )/V2 3 

Insist A^ = (A%^-A a L)i )/V2 3 

V^V Xf =X a ? +X^ a ' 3 3 -3 

iv a l5 ^ X°f = i(X a f i-X^P) 3 -3 

^fa^ Bf^ = (B^+B^)/V2 3 1 

^ B% v = i [W* - Blf) /V2 3 1 

Although the flavor representation is correct, in four dimensions a complex two-form has too many degrees of freedom, 
the reason is that the tensor operators are not all independent, but satisfy the duality condition (JTJ. This relation 
implies that the complex two-form has to be imaginary anti self-dual. 

Defining Fl and Fr as the field strengths of Al and Ar, H = dB — iA^ A B + IB A Ar as the three-form field 
strength of B, and DX — dX — iA^X + iXAr as the covariant derivative of X, the action takes the form 

S = J d 5 x^tr[- ^(F 2 + F 2 ) + g 2 x (\DX\ 2 + 3\X\ 2 ) - 2g B (i~(B A fft _ B t A H ) + m B \B\ 2 

+ ^(X^FlB + BF R X^ + h.c.) . (23) 

The trace is taken over the gauge indices. The factors of the AdS radius I have been absorbed in the coupling 
constants or the masses. 

Let us comment on the different terms. The first term is the kinetic action of the gauge fields, its coefficient was 
fixed in the original hard wall model comparing the expansion of the holographic vector-vector correlation function 
at large Euclidean momentum with the OPE of QCD The result was 

- = -^2- (24) 
gt ~ 12n 2 ■ {M) 

The second term is the scalar action, that is usually canonically normalized, which can be achieved by rescaling 
gxX — > X in lf23|) . The asymptotic value of the scalar field close to the boundary z — determines the quark mass 
m and the condensate (ipi/i) in the dual theory. With canonical normalization, gx now appears in this relation: 

= ^V(*)12X2. (25) 




The reason that there are not two independent normalization constants in this relation is that the expectation value 
of the mass operator (tpip) can be obtained from varying the on-shell action with respect to m. The value of gx 
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was determined in ref. [7j comparing the holographic scalar-scalar correlation function at large Euclidean momentum 
with the OPE of QCD, 

5 | = gr. (26) 

The third term is the action of the two-form held. The kinetic term in the action has been replaced by a Chern- 
Simons term, and the mass niB is a free parameter. We will see latter that the right self-duality condition for the 
two-form field can be derived from the equations of motion of this action by fixing the mass. Then, following the 
usual procedure, we will compute the holographic tensor-tensor correlator, expand it at large Euclidean momentum, 
and match with the OPE of QCD. The last term is the most general gauge-invariant term of dimension five or less 
that couple isospin triplet vector mesons and preserve parity and charge conjugation in the dual theory [9J. 
We can rewite the interaction term in ([23)) using real fields 



tr 



(X*F L B + BF R X^ + h.c.) = itr (x + ({F v , B+} + i[F A , + A_ ({F v , B_} - i[F A , £+])) . (27) 



The term (|27p determines how chiral symmetry breaking affects to the isospin triplet 1 vector mesons. Without 
this coupling the spectrum will be determined by the equations of motion of the V field, but once we introduce it, 
the B + field and the V field are coupled. 

Using the expression (1231) for the background scalar field and taking the trace in the action ([2"3")l . we get in the 
vector sector 

S V = fd 5 xyT9[--k E FiMNF iMN + ^e MNLp Q{B- MN H +LPQ - B +MN H- LPQ )- 

-gems J2 B «mnB% in + ^v(z)F v mnB^ n ], (28) 

£*=+,- 

where we have suppressed gauge indices. In total we have introduced three new parameters, gs, ms and A. We will 
now fix ijig imposing the self-duality condition and gs and A using the matching with the OPE of QCD. The only 
free parameters left in the model are the mass m, the condensate a = (tpip) I Ox ano - the ^ scale 1/ z m . Although 
by introducing the B field we have added a new sector of vector mesons and therefore of masses and decay constants 
we can compare the model with, we have not increased the number of free parameters. Notice that the axial sector, 
involving the fields X and A is untouched. 



IV. EQUATIONS OF MOTION 

Our next step is to calculate the equations of motion for the fields B± ,B±* , and V M from the action (f28|) . We will 
write explicitly all the factors involving the radial coordinate and raise and lower indices with the flat metric, using 
gMN = jjVmn- Greek letters for the indices will refer to the flat Minkowsky directions and capitalized italic letters 
will include the radial direction z. Let us consider first the case with no interaction, A = 0. From (f2"5| we get the 
equations of motion for the components of the two-form 

± \e MNLPQ H TMNL + m B B PQ = 0. (29) 

Including explicit powers of z we have: 

±e «*f>*» H + = 0, (30) 

z 

± e ™0yn H + *m}B B g = 0; (3i) 

z 

with the definition e za p^ v = tafiiw =>■ e zal3f2l/ = —e al3tJ " L ' . After some algebra, one can write the equations of motion 
as 

zd z {zH^ u ) + z 2 8 a H^ + = 0, (32) 



■2 



Az 2 



B»± = 0. (33) 
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Note that equation (|29j) is of the type dB — *B = which implies the constraint d * B — 0. In components, 

dpBg = 0, (34) 
d„B^ + zd z ^B z ± " = 0. (35) 

Using (1351) in (|32|) and (|33[) we can eliminate d z B± z from the former and B± from the latter. Expanding solutions 
in Fourier modes of four-momentum we get: 

7?dlB^ + zd g Bg + (z 2 q 2 - ^ = 2izq^Bf, (36) 

d 2 z B£ z - -d z B% + (q 2 + B£ z = 0. (37) 

Let us use the following decomposition 

B£ = iq^T^} + i€^ a Pq a T ± p . (38) 
If the two- form field is dual to an operator of conformal dimensions A = 3, then its expansion at small z should be 

T£ = -T^l + z log zT^ ± + zT^ ± + ■■■ 

T± = -T ± + zlogzT ± +zT ± H (39) 

This is possible if we set m 2 B = 4. Setting = 2 and using the original equations ([30]) and (|3T|). from the leading 
~ l/z 2 term we get the conditions 



? (o)^ = ±T (o)^ (4Q) 



The conditions (14TJ1) above imply that 



B^ + , v + \e, va[i B^ = . (41) 



Let us define fe M „ = B^ + ^ v — iB^> In terms of b, the condition above means it is imaginary anti-self-dual 

6 + i*6 = 0. (42) 

The conjugate w is imaginary self-dual. 

Let us now study the effect of adding to the action a kinetic term for the two-form field, with a relative coefficient 
C> 0. 



AS = g B C fd 5 x^[ J2 H aM NLH^ INL ], (43) 
«=+,- 



Given arbitrary mass ms, we find the following equations 

CV L H ± LMN ± J £ MiVABC J ff T Asc + m flJ B ± MiV = . (44) 
Where Vl = -j=dLy/—g- Assuming that the leading term in the small- z expansion of the two-form field is B±" ~ 



z v B^^± , one finds that B^ 1 ^ satisfies an imaginary (anti) self-duality condition if 

Cv 2 T2v-m B =0. (45) 

The upper sign corresponds to the condition (|4"Tj) . The derivation is valid if raj — CV 2 ^ 0, otherwise the leading 
term is a logarithm and the analysis is different. The equations for plus and minus components can be decoupled, 
giving a fourth order equation 

(4 + 2m B C)V K H ± Kp Q + m 2 B B ± PQ + ^-V L [{g KP g c ^g DL + (PQL)) d K (gcufjDvV S H ± SUV )] = 0. (46) 
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Where by (PQL) we denote permutations with a minus sign if they are odd with respect to the first term. To leading 
order in z, the equation for the spacetime components B± v imposes a constraint on v in the form of a quartic equation 



CV- {A + 2m B C)v 



= 0. 



(47) 



We can rewrite this equation as 



(Cu 2 + 2i>- m B )(Cu 2 -2v- m B ) = 0. 



(48) 



Therefore, when the solution is (anti) self-dual the quartic equation is automatically satisfied. This shows that we 
can always impose the right self-duality condition, although we can also have solutions with opposite self-duality 
conditions if we change we have to set those to zero by hand. Instead, we will drop the kinetic term, we will 
show that it does not affect to properties of the model like the meson spectrum, although it can affect to correlation 
functions because it contributes to the boundary action. 
If we set C = 0, then the equations (|16")) become 



KPQ 



-B 



PQ _ 



= 0. 



(49) 



Let us now do the following trick, we can rewrite (|44[) as 

1 



C 



LAIN 



-B 



AIN 



AIN ABC 



H 



tabc 



m B 



m 2 B C 



fiMN 



(50) 



Then, if 



m B = m B 



(51) 



the solutions to the C = decoupled equation (149j) with mass rh B would make the term that multiplies C in the 
first bracket vanish, so one would recover the C = equations again. Now let us fix the asymptotic expansion (|39[) 
(v = —1) and impose the self-duality condition (|41l) . By fixing the mass to the value m B = 2 + C we can solve the 
system with C ^ using the solutions for C = 0. 25] If one examines the equations involving the interaction term 
below one sees that the same is true if the coupling is rescaled appropriately. An exception to this rule may be the 
case m B = 0, where there is an additional gauge invariance associated to the two- form field SBmn = d^M-^-N] and the 
separation in two parts of the equations of motion involves introducing gauge non-invariant terms. It also coincides 
with the case v = 0, where the leading solution is logarithmic. We will neglect this case and set C = from now on. 

We now consider the interaction term, it does not affect to the leading asymptotic behavior, so the value oim B =2 
is not changed. The equations of motion of the B + two-form are 



z 80s z 



d a zH_ 



avz , B? A v[z) 



zd a H^ v + d z zH z _^ 

d a zH°L vz + — = 0. 
z 



z 



13 Z 

B» v 



V 1 



r v 



8g B 



v'{z)F VaP e a ^», 



(52) 
(53) 
(54) 
(55) 



The equations of motion of the vector fields and the constraints for the two-form fields are 



d z -B v + z 
z 



z 



d z -B v J 

z 



d B^ 



z 



A 
8^ 

A35 



= 0, 



z 

d z V -^Bl*- 



22 v2 ( z ) 



v(z) 



z 

v(z) 



d F^ 



8 B^ v 



9 5 9 x 



(56) 
(57) 
(58) 
(59) 
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We now expand in Fourier modes and simplify the system to: 

d z B v + z - 44 B + - iQ» B T = 



2-iru 



div 



f(z) 



A v'(z) 

^9b g(z) 



d z V v 



9(z) 



d z v v + q Z V u = 



g{z) 

P a A" - -d z A^ + q*A» = g 2 5 g 2 x V -^^ 



d z B v _ z - -B U J - iq u B^ = 

z 



(60) 
(61) 
(62) 
(63) 



Where f(z) = - + \ zv(z)(^-p-)' , g(z) = 1 — x v ( z ) 2 an d X = ^ 2 9 2 /{%9b)- Note the structure of the equations, we can 
understand how the fields mix among themselves by considering parity conservation: the vector mode V^ L (negative 
parity) couples to its vector partner B+ z (|6"Tj) and to the tensor component of negative parity (I5"2"j) . We can further 
decouple the system of equations above: 



d z - (C\(z) - q 2 ) 



B V J = 0, 



where 



B1 



9(z) 



C x {z) = d.. 



m , 1 , Xv' 2 
g(z) z 2 g{z) 2 '' 



9(z) 

The equations for the tensor components of the two-form field are 



9 2 (z) 



(64) 
(65) 

(66) 
(67) 



[zd z zd z - 1 + z 2 q 2 } B^ 



A 



} ( z )ql^ - ^ld z q^ 
9{z) 



+ z 



[zd z zd z - 1 + z 2 q 2 } B^ = - l ^zv'(z)q [a V p] e a ^ u + 2ziq^B v } z . 



9(z) 



: j iq^B^ z , (68) 
(69) 



Our final equations of motion can then be divided in the vector (|61| . axial vector (|63[) . and two- form components (|64l) . 
([65]). (|68| and (69]). A more convenient grouping is in four decoupled sets: {A^}, {V, B^ z , (5%-P£)B% u , P£PgB°f}, 

{P$P$B ( f} and {B^-, (6%- P£)B°:»}. Normalizable solut ions of the first two sets correspond to 1 ++ and 1 mesons 
in the dual theory, respectively. The last two are not independent since they are coupled in the original system of 
first order equations (|30l) . (I5T|) . and normalizable solutions correspond to 1 H mesons. Notice that b^" = zPj^VpB^f 
satisfies the same equation as B^_ z , and that this one is the same as the equation for vector mesons (|6ip in the absence 
of the interaction term A = 0. Therefore, the interaction lifts the degeneracy between 1" and l" 1 mesons. 



A. Boundary expansion 



We now proceed to do a Frobenius expansion of solutions close to the boundary at z = 0. This will be useful for 
both the calculation of renormalized two-point functions and the calculation of the meson spectrum. Using (1381) in 
(1501) . we find the conditions 



T(d z T T ^ + B 



d z T^_ + ~T*t = 
z 



z 

9b z 



1- 



0. 



d z Tl + -T» = 0. 



(70) 



10 



And from (13TT) we have 



-q^_ + -B^ = -^-d z V^ 

q 2 Tl + -B» z = 0. (71) 
The expansion of the vector components at small z is given by (|39[) and 

A^ ± = A^l + z 2 log zA^l+z 2 A^l + - ■■ , 

Bg = B^ ± + z 2 \ogzB^l + zW ± + ■ ■ ■ , (72) 

where we have defined A+ — V 1 and At — A^. Expanding ([70]) . (|7T|) . (|57|) and ([55)1 for small z we find a set of 
conditions that allows us to solve for the coefficients of the logarithmic terms and give us a relation between the 
leading terms, dual to sources in the field theory. In particular we recover the imaginary self-duality condition (|40j) 
for the components of the two-form field. Defining A+ = A/(8<?_b), A_ = 0, q\ — q 2 and q 2 _ = q 2 — g\g 2 m 2 , we can 
write them in a compact form 

B<®1 = ±q 2 T^ ± , T {0) " T =±T^ ± , 
T (i)P = ^ T ( 0) M _ h± mA ^ T (1) T = T |t(») M ± ± ^mA^l, 

T% = T T(% ± \{q 2 T^ ± + ~X ± mA^l) A^ ± = -\{q 2 ± A^ ± - g 2 X ± mq 2 T^ ± ). (73) 

V. HOLOGRAPHIC RENORMALIZATION 

We will follow the usual holographic procedure to compute correlation functions, deriving the on-shell action with 
respect to the sources of dual operators. The action usually diverges, so we will introduce a cutoff at a small value of 
the radial coordinate z — e to regularize it. We will introduce counterterms following the usual prescription [16[ to 
make the action finite before removing the cutoff e — > 0. 

The on-shell regularized action is 



S„... = J d A x^~g ^sTgT E F±'<*A±v- WaT B + , Z »A + , V 



«— j z =e 



Scs- (74) 



Where we have introduced the cutoff e and the overall sign corresponds to taking the lower limit in the z integral. 
Scs is the contribution of the two-form Chern-Simons action. The action has the form 

Scs = ~9b J d 5 x l - tr (B A £T f - A H) = 9 -j / d 5 xe MNLPQ (B_ MN H +LPQ - B +MN H_ LPQ ) . (75) 

A variation gives, to leading order, 

SScs = - y I d 4 x (B-^SB+np - B + ^5B_ afi ) . (76) 

Where the fields are evaluated at the cutoff z — e. The condition (|4Tj) implies that we cannot vary B+ and B- 
independently, but since B+ + *_B_ = at the boundary, we should treat B+ — *B- as the variable boundary value. 
In order to have a consistent variational principle we need to add a boundary term to the action, of the form 



S = — 



J d 4 x ^tr (r^Bl^) , (77) 

where the indices are raised with the induced boundary metric 7 MZ , = s~ 2 rj^ v . The variation of this term, with expicit 
e factors, is 

SSq = ^J di < B r sB +^ + BJ^SB-^) . (78) 
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The sum of the two variations gives 

6(S CS + S o) = ^j d 4 x(B + ^ + l -e^B_ a0 )8(B +llv - \e^B. 



up 



(79) 



This gives a consistent variational principle. 
The total on-shell regularized action is 



So.,. = [d 4 x -~]T d*KA a » + X-B^V, + 9 -^J2 B *^ B « V 

Z9 5 Z a=± Z 6 a=± 



(80) 



The bulk contribution vanishes on-shell. 

Expanding for small e we find that, for Fourier modes of momentum q^, 



-d T A1A 



a=± 



(2A^ ± , + A^ ± ,)A^ ± + 
+2log(Qe)A^ ± ,A^ ± , 

Vlog(Q £ ) [q 2 T^ + T^ + , + q 2 ¥ 0) l¥ 0) + , _^ mT (0)^(0) + ^ + 

+ 8<z 2 (r (0)^ T (2) + Aj _ T (0)^( 2 ) + Aj) + ... 



(81) 
(82) 



(83) 



Where the dots refer to local terms in the sources, they will not be relevant because we can remove them with finite 
counterterms. 

To this action we have to add some boundary counterterms to remove the divergences that appear as e — > 0. 
As expected, the leading divergence 1/e 2 docs not appear in the action of the two- form field. There are however 
additional logarithmic divergences. In order to completely cancel them we need more counterterms, of the form H 2 , 
F 2 , (dX) 2 and XFB. More explicitly, we have that the finite regularized action is 



Sreg — S Q . S + S\ + D2 + S3 _|_ + S3 



1S4 + Sfinitc 



where 



Sl=Cl [ d 4 xl0g(»s)^Yl H -^ H a V ° ', 
J a=± 

S 2 = c 2 J d 4 xlog(ps)^X + F v ^B^ , 
S 3 ,±=c 3t± J d i x\og(ne)^F±^F£ v . 
S 4 = c 4 J d A x\og{^e)^{D^X)\D^X). 



(84) 

(85) 

(86) 
(87) 



Notice that we can also have finite conterterms, corresponding to Si, S2, S3 ± and S4 with no log factors. We will 
introduce the finite counterterms in Semite and use them later on. 
In the e — > limit, the counterterms become 



51 r 

5 2 ~ -2mq 2 A^ + fl T% 
S 3± ~ -2q 2 A^ ± ,A^l 

S 4 ~ m 2 A^_A^_^. 

One can cancel the quadratic and logarithmic divergences if 



6(q 2 ) 2 (T^ + ,T^ + +¥°\,T W +), 



(89) 

(90) 
(91) 
(92) 
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Now one can take the e — > limit, and use finite counterterms to remove the pieces that are local in the sources 
d 4 q 



S r 



(2nY 



+ ^£,2(^(0)^(0)^+^(0)^(0) 



(94) 



A. Two-point functions and matching to QCD 

We now proceed to write the renormalized action in terms of general sources w M and t^. The transverse vector 
field is 



A<?>\ = P/*a c 



m lilt , 



The transverse tensor is 



The longitudinal part of the tensor is 

2 

Using the expansion 



we indeed find 



t 



1 



_ p " p Pj. _ 

M ?y ~~ 2 ^ "1 ^ ~~ 2q 2 Ml/ 



1 



t 1 - if a a T""' t u - ia< T (0) 1 
We now assume that the subleading terms are proportional to the sources 

A^ + , = G vv { q 2 )A^ + , + G VT { q 2 )T^\^ 
T^ + , = G TT { q 2 )T^ + , + G TV { q 2 )A^ + „ 



7f&) , _ n + ,2^(0) 



T 



— G TT ( q )T 



And using ([T 



^(0)^(0)A 



? (0) ? (0)M 



1 ~ - r fl.r U 1 ■ jpMU;aP. 

'4(^)2^ 



r9 r a/x9 r ^ 



a/3 



(g 2 ) 2 



4(<Z 2 ) 



2g 2 

9f-iaQisQ/3 t. apt y§ 



4(f 



(95) 

(96) 

(97) 

(98) 
(99) 



(100) 



Introducing these expressions in the renormalized action (|94p and deriving twice with respect to the sources, we find 
the following correlation functions 



iCv a V) = s ab W - ^ F )n vv (, 2 ), 



n 



TV 



VT 



( Q 2 ) = id^Wqr-^PqnUvTiq 2 ), 

i/.U- 



n^ Q/3, °V) - s ab n+ T ( q 2 )F^' al3 + s ab n^ T ( q 2 )F^' ali 



(101) 

(102) 
(103) 
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Where 



nyy(q 2 ) = - 


l Q 2 


2 2 
+ 2 „GVv(<7 ), 


(104) 


nvT(<7 2 ) = -^mloj 


Q 2 l 
5 M 2 2 5s V 


4 

G V T(q 2 ) + ^dBGrviq 2 ), 


(105) 


n+ T ( g 2 ) = 


3B, Q 2 

-T l0 V 




(106) 



n TT (g ) = -- -log— -- — G TT (q ) 

Comparing with the expressions (p~8|) . (fl~9|) and (|20|) we get 

1_ _A^_ JV^ 3A e 

.gf - 12tt 2 ' 5S " 16^ 2 ' 4^ 2 ' 

Together with (|26|) . this fixes all the coupling constants of the bulk action. 



(107) 
(108) 



VI. MATCHING TO THE MASSLESS THEORY 



We have used the coefficients of the logarithmic divergences of the correlation functions to fix the parameters of 
the model. Notice that in the massless limit m —> the logarithmic contribution to the vector-tensor correlator 
(|105p vanish. In QCD, perturbative contributions vanish to all orders, so the only contributions left come from 
non-perturbative physics, this is clear in (|20j) . where the leading term when the mass is zero is proportional to the 
condensate. Therefore, for massless QCD this is the term we have to match to fix the value of the parameter A. Since 
the coefficient of this term is independent of the mass we should get the same value for A, we will se that this is indeed 
the case, so the model passes this non-trivial consistency check. 

In order to find the non-perturbative contributions to the OPE, we need to compute the functions l G(q 2 y that 
appear in (|100[) and plug them in the expressions for the correlators that we have found in the previous section. 
The overall strategy will be to solve the relevant equations of motion and match the near boundary expansion of the 
solutions to the coefficients of the series defined in (|39j) and (|72|) . 

As we mentioned before, there are four coupled equations describing the negative parity mesons. Plugging (|38j) in 
equations (f60|). (j6l"j) , (l64|) , (|68|) and (l69l) and keeping only the negative parity modes (V,T",T_,B^"), we have the 
equations: 



^-M dz -{C 1 (z)-q 2 ) 



9(z) 



Bl z = - 



&g~E 



d 2yu /(*) 



9(z) 



d z V v + q*V v = ^MHb v 



+ ' 



c a (*)v* + ^W-v 



[zd z zd z -l + z 2 q 2 ]T» = - — 

8g B 

[zd z zd z - 1 + z 2 q 2 ] T_ = --^zv'(z)V v 



9{z) 



9(z) 



! j z 2 44+ :)B?, 



(109) 
(110) 
(111) 
(112) 



With the aditional constraint: 



f(z) 

d z B\* - J -±iB v + z 



iq 2 Tl 



A v'(z) 

8.9s g(z) 



d 7 V v . 



(113) 



The mixing term proportional to \v(z) in the equations of motion is a small perturbation for small values of z, as 
v(z) falls off towards the boundary. As the large Q behavior of spatial correlators with Q 2 — —q 2 is dominated by the 
small z behavior of the solution, we can determine the short distance behavior of correlation functions analytically by 
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treating A as a small parameter and solving the the equations of motion perturbatively in A. However, we don't need 
to solve all the four equations. First, the constraint above implies that B vz is not independent, and the relations ([75)1 
imply that close to the boundary, TY is not independent of T_. For convenience, we will focus on equations (|111[) 
and (|112p . Dropping all terms of order A 2 in the equations above, and taking the ansatz: 



V = V (xX + Vx(x)W 
T v _=T {x)t v _+T x {x)v v 



(114) 
(115) 



Our problem is reduced to solving the equations below 



d*--d x + <f)V (x)=0, 

x 



di - -d x + l J V x {x) = A(ai + a 2 x z )B {x), 
[x 2 d 2 x + d x - 1 - x 2 } T (x) = 0, 
xdl + 8 X - ±±E- T\(x) = -X(T lX + T 2 x 3 )Vo(x). 



(116) 

(117) 
(118) 
(119) 



Where ai 



-hr, OL 2 



34 7 1 1 



and To 



3<T 



Vq, Tq are the homogeneous solutions and V\, T\ the 



Q* > ■■•1 — 8g B Q — 8g B Q 3 ' 

perturbative corrections of order A. Bq satisfies the same equation as Vq. Focusing first on the vector mode, Vo and 
Bq have well known solutions in terms of Bessel functions: 



V Q (x)=x{ah(x) + K 1 (x)), (120) 
B {x)=x{bh{x) + K 1 (x)), (121) 

Where a and b are constants fixed after imposing the IR boundary condition. For Vq(x), following previous 
work we are going to choose a Neumann boundary condition at x m , d x Vo(x)\ x = which allows us to set 

a = ^"(x™) • ^° cnoose the appropriate boundary condition for Bq we note that a Dirichlct boundary condition 

imposed on Bo(x), which sets b = — ^^ x m -j , implies, by use of the constrain (|113p to leading order, a Neumann 
boundary condition for the leading tensor mode, but a choice of Neumann boundary condition for the former is 
not consistent. As Bq and Vq describe completely independent modes, this already captures the most general solution. 

To compute V\ we will use a Green's function method. Note that we can write a solution for equation ()117[) of the 
form: 

V\{x) = A f ™ dx'[ ai + a 2 x' 2 }Bo(x') Gv{x ; X '\ (122) 

J x e x 

Provided Gy satisfies the equation: 

d 2 x - ^d x - lj G v (x,x') = xS{x-x'). (123) 

With boundary conditions Gv(x e ,x') — G' v {x,x m ) — 0. We solve the equation above in the two regions x > x' and 
x < x' and match the two solutions at x = x' . It is not hard to show that Gv(x,' x) can be written as: 

G v (x,x') = AD _ BC [Ah(x > ) + BK 1 (x > )][Ch(x < )+DK 1 (x < )}. (124) 

Where .t< > — {ruin, max} (x, x') is book keeping notation to specify the two branches of the Green's function. The 
coefficients are A — — K a (x m ); B — I Q (x m );C = K 1 (x e );D = —Ii(x e ). Taking the limit x e — > we can set above 
D = and C = 1. Replacing back in (|122l) we have: 

rx 2 rx m 

V\(x) = -A / dx'[aix' + a 2 x' 3 }K 1 (x')h(x') - A— / dx'[a x x' + a 2 x' 3 ]Kf(x'). (125) 

2 J x 
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Where we used: 

x{Ch{x) + DK 1 (x)) 
AD-BC 

x(Ah(x) +BK 1 (x)) 
AD-BC 

and 



= -xh(x) ~ -— , (126) 
~ -1, (127) 



Bo ~ xifi(x). (128) 

We also have ignored all the terms proportional to -g ~ e -2a m smce m the limit of large momentum these terms 
vanish quickly. Physically this means that as the momentum increases what happens in the IR region becomes less 
important, as expected. In fact, these solutions near the boundary and for large momentum become oblivious of 
the IR boundary conditions, since both Neumann and Dirichlet boundary conditions will enforce factors that fall off 
exponentially. Moroever, we will ignore the contribution of the first integral, that is negligible since these contributions 
vanish too quickly near the boundary. 



We then have: 



y A (x) = -Ay^ m dx'(a lX ' + a 2 x' 3 )K 1 (x')K 1 (x , )~ Aa X y Q+logz) -Aya 2 . (129) 



The near boundary solution for V u is: 



V v (x) = V (x)v v + V x b v = [ 1- — + — logx)v v + ( Aai— ( - + log a: I - \—a 2 ) b". (130) 



Ag = V™ = -^y(o) - AQ 4 - f ) Tf . (131) 



Matching the solution above to the expansion defined in (f72j) . and using (|73|) . we find: 

4~' V4 3 

Therefore: 

G + VV {Q 2 ) = -^, G+ T (Q 2 ) = -AO 4 (f-f)' (132) 

Following similar steps, we can now compute near boundary solutions for equations (|118l) and (|119[) . The homogeneous 
equation has a well known solution of the form: 

T-{x) = ch{x) + Kx{x) ~Kt(x). (133) 

For the second equation, again, we can write a solution with a Green's function: 

T x (x) = -x[ dx'{T lX ' + T 2 x' 3 )K 1 (x')G T {x,x'), (134) 



where analogously to the previous calculation, Gt(x,x') satisfies: 



xdl + d x 



l + x 2 ' 



G T (x,x') = 6(x-x'). (135) 



It can be shown that: 

G T (x,x')^-^i(a;>)/i(a;<), (136) 
in the limit where x e — > and Finally, we find that the solution for the tensor field is: 

Tx(x) = A| £ dx'[{T lX ' + TW^Kxix')}^') =\^l-^ x , (137) 

therefore, 

T (2), = _Q^ T (o), + Q m TA V(0K (i3g) 
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However, we are really after T( 2 ^ = G T t(Q 2 )T {0 ^ + G^ v (Q 2 )V^^. To compute the latter, we use the relations we 
have found previously: 



Solving for t} 2)ai we get: 



-.(2)/i 



-T 



_A_ 



So the result is, 



Q 2 
4 



Gj T (Q 2 ) — —-7-) G^ y (Q 2 ) 



5Ar 



5Ar 



32G B 32.g B Q 2 



y(0)M 



(139) 
(140) 

(141) 
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We can now compare with the OPE of vector and tensor correlators 
nonzero contributions are 




G 



Q 2 

vv = ~~r> Gvt = Agfa, 

r-vk _ Q 



TV 



Xa 



>Q 2 



(142) 

Setting the mass to zero, the only 

(143) 
(144) 



The contributions Gyv and Gtt give contact terms that can be removed using counterterms in the regularized action. 
The only non-pcrturbativc contributions to this order are 



1 2\ Act 



6g% Q 2 



(145) 



Comparing with the OPE in QCD (|20|) . we see this term has the right coefficient, including the sign. The results we 
have obtained are insensitive to the details of the IR, so they should be valid for any models that are asymptotically 
AdS space. However, the value of the condensate itself and other quantities like the meson spectrum will be sensitive 
to IR physics. In the next section we will study how the inclusion of the new terms in the action affect to some of 
these quantities. 



VII. MESON SPECTRUM 



So far we have discussed the UV physics of our model, focusing in the matching with the OPE of correlators in 
QCD. We will now comment on some of the IR physics, in particular the meson spectrum. In our analysis we have 
seen that the two-form field splits in a transverse part and a longitudinal part, that mixes with the vector fields. We 
can summarize the correspondence between the fields and meson states in the following table: 

B^ v mixes J PC mesons 
transverse — l" 1 
longitudinal 1 

The lightest isospin triplet states that can be found in the Particle Data Group (PDG) review [13], are 

meson J PC mass (MeV) 

6i(1235) 1+- - 1229.5 ±3.2 

,0(770) 1— - 775.49 ±0.34 

<zi(1260) 1++ ~ 1230 ±40 

Notice that according to the PDG estimate, the 61 and a\ mesons are almost degenerate, although the error in the 
estimate of the a\ mass is very large. Other estimates give a mass to the a\ ~ 1255 MeV, with somewhat smaller 
errors [l8l |. 
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In the holographic model the b\ state is obtained from the transverse components of the B field, that are decoupled 
from the rest of the fields. The degeneracy between the p and the b\ is broken in our model, thanks to the interaction 
term proportional to A in ref. (|23|) . Had we not considered this term, the spectrum would be degenerate, as has 
been observed in ref. Q. So we should include this cubic interaction term both from the perspective of the large 
momentum OPE and from the properties of the meson spectrum. 

We follow a similar procedure as in ref. [1] to compute numerically the lowest masses of the vector meson spectrum. 
We must specify suitable boundary conditions for the fields at the IR radial cutoff z = z m (Neumann or Dirichlet) 
and at the boundary z = (normalizability). Solutions do not exist for any value of the four- momentum q 2 , but only 
for a discrete set of values, which correspond to the masses of mesons in the holographic dual m 2 = q 2 . We have 
checked that our results for the meson spectrum and the pion decay constant /„• coincide with those of ref. [1] when 
we set the coupling A = 0. 

We start with the spectrum of 1+" mesons, dual to the field components {P£P£B°f} and {B^ z , (6£ - P^B™}. 
Notice that we can solve first for B^ z in (|65|) and then use (|69|) to solve for (8£ — P£)B" l/ . As we have explained (|68|) 
is equivalent to (p5|) , so for the purpose of finding the masses it is enough to focus on (j6"51 . Close to the boundary, a 
normalizable solution has the asymptotic expansion (j?2")) with B^ 1 ^ = 0. At the cutoff we impose Neumann boundary 
conditions, since for Dirichlet boundary conditions there is a normalizable solution at q 2 = 0, which would be dual to 
a massless vector meson. Normalizable solutions are Bessel functions Bt z = b^ zJ\{\q\z) and the Neumann boundary 
condition is satisfied for values of the momentum such that Jo(\q\z rn ) = 0. Then, the mass of the lowest mode is 

m bl z m ~ 2.405. (146) 

Notice that this value is independent of the quark mass and condensate. The remaining modes do depend on them 
and we have to solve numerically the equations. 

We will first solve for modes dual to pseudoscalar mesons, whose lowest mode corresponds to the pion. We need to 
solve the set of equations [l[ 

v"-~v' + 9hx^(*-v) -0, (147) 

_ g y + ^4^' = 0. (148) 

z z 

For this, we first derive a single second order equation by solving algebraically for 7r in the first equation, plugging the 
result in the second equation and defining (ft — (p 1 . Then, using g\g 2 — 3 and defining h(z) — 3v(z) 2 /z 3 , we obtain 

<ft"{z) + j^(z) -{['-) - , - + zh-q 2 } 0(z) = 0. (149) 

Normalizable solutions at the boundary behave as <fi(z) ~ z and we impose a Dirichlet boundary condition at the 
cutoff for the field (ft. Then, for given values of mz m and az^ we find the lowest value of qfz^ = m 2 z^ n such that a 
solution satisfying the boundary conditions exists. We can then use the physical value of the pion mass ra, = 139.6 
MeV to fix the scale z m . 

The spectrum of axial vector mesons 1 ++ can be found by solving equation (|63|) . From (f72j) a normalizable solution 
^4(0)^ _ q ya-nigheg a ^ ^ ne boundary. At the cutoff, we impose Neumann boundary conditions. Finally, the spectrum 
of vector mesons can be computed from the system of coupled equations (pTlj) and (f6"4")l , with conditions A^ 1 ^ = 0, 
— in the expansions at the boundary (|72|) . Regarding the boundary conditions at the cutoff, we must be careful 
since equations (|64|) and (|68|) have an additional singular point at z* such that g(z*) — 1. We are then constrained 
to values of the quark mass and the condensate such that z, > 1 or to impose suitable boundary conditions at the 
singular point. A quick analysis shows that the two possible behaviors of solutions close to the singular point are 
~ (z* — z) 1 ! 2 and ~ 1 for and ~ [z — z*)( 1±v/ ^-'/ 4 for B^ . We can then make the solution regular by imposing a 
Neumann boundary condition for and a Dirichlet boundary condition for B^ z . 

For the values of the mass and the condensate we have explored 0.0001 < mz m < 0.1, 0.0125 < az^ < 0.5 we do 
not find a realistic spectrum of mesons, the lightest vector meson 1 is always heavier than both parity even mesons 
1 ++ , l + ~. For larger values of the mass we can understand this as a consequence of the singularity at z = z*. The 
curvature of AdS makes the classical problem of finding normalizable modes effectively as the quantum mechanical 
problem of finding the energy spectrum of a particle in a box, with one of the walls at the cutoff. For the 1 modes 
we are forced to impose boundary conditions at the singularity z» < z m , so the "box" is smaller and the spectrum 
is lifted to higher values. This could be a problem of how infrared effects are implemented in this particular model, 
maybe different constructions like the soft wall could avoid this issue. 




18 



There is a way to find a more realistic meson spectrum, with the parity odd vector meson below the other modes. 
Instead of introducing the cutoff, we can impose boundary conditions for the vector mesons at the singularity even 
when it sits at a radial position beyond the cutoff z* > z m . For large enough values, the vector mesons become 
lighter and the spectrum can be tuned to realistic values, for instance for mz m = 0.0005, crz^ = 0.05375 we find 
that m p ~ 753.95 MeV, m at ^ 1238.24 MeV and m bl ~ 1237.87 MeV. Although this would fix the meson spectrum, 
there are other quantities that are important to determine whether the model is phenomenologically viable. One such 
quantity is the pion decay constant, f„, that in QCD is approximately / w ~ 91.92 MeV. In the holographic model it 
is given by the formula [1] 

2 _ 1 d z A(z) 

where A(z) is a solution to (l63l) satisfying A(e) = 1, A'(z m ) = 0. With the parameters that give a realistic meson 
spectrum, the value of the pion decay constant is quite low / w ~ 4.07 MeV. 



(150) 

Z—F. 



VIII. CONCLUSIONS 



We have carried out, for the first time, a complete treatment of the hard-wall model including all fields dual to 
operators of free field theory dimension 3. We followed the standard procedure of fixing bulk parameters by matching 
the short distance behavior of correlation functions to perturbative QCD. Reassuringly, the structure of the correlators 
we obtained from our holographic model precisely matched the expressions in perturbative QCD, so this program 
can be carried out consistently. With this matching in hand, we calculated physical properties of mesons which, 
unfortunately, no longer match QCD. While this result casts into doubt whether the simple hard wall model can serve 
as a good stand-in for QCD, one may hope that an improved IR model could potentially lead to a better spectrum. As 
our analysis of the short-distance behavior of correlation functions only relies on the UV asymptotics of the geometry, 
the action we derived (including the numerical values of the coupling constants) should serve as the starting point 
for any such exploration of complete (in the sense of including all dimension 3 field theory operators) holographic 
bottom-up models with alternative IR boundary conditions. As we discussed in section IIV1 it is possible to modify 
the bulk action of the two-form field by adding a kinetic term, giving a one-parameter family of theories with the 
desired self-duality condition and asymptotic behavior. Since this will modify the boundary action, in principle the 
value of the bulk couplings will be shifted when the matching to QCD is done. It is possible then, that by changing 
this parameter, a more realistic spectrum can be found. 

Let us point out some differences between our approach and what one expects in a top-down models like Sakai- 
Sugimoto [H, [2(j , based on a string theory construction. The matter content of the model is such that it coincides 
with large-A c QCD at low energies in some region of parameter space where the UV theory is weakly coupled. In 
particular, 1^ mesons should be part of the spectrum. However, in the holographic description where supergravity 
is valid such modes are missing. This should not come as a surprise: since the tensor operator is not a BPS protected 
operator, its conformal dimension can receive large corrections of order ~ A 1 / 4 , where A is the 't Hooft coupling. In 
the holographic description this means that the tensor operator is dual to a field with a mass of order of the string 
scale, and therefore beyond the supergravity approximation. Since corrections to non-BPS operators are very large, it 
is even possible that the lowest l" 1 meson is not described by a field dual to the tensor operator we have considered 
in our model, but to a different operator with the same quantum numbers but larger conformal dimension in the free 
theory. This indeed seems to be the case in the Sakai-Sugimoto model, where the l" 1 mode is described by some 



components of a symmetric field in the bulk [21( . Clearly, in this case we do not expect that the OPE of the model 



will match with that of QCD, so in some sense the approach of refs. HQ is closer to the top-down model. However, 
if the dimension of the tensor operator is chosen to be larger than 3, it is more difficult to argue that the effective 
theory description in the bulk stays valid anymore. 

We have studied the extension of the model that takes into account l" 1 mesons, like b and w. In principle the 
model can be further extended to include other modes in the QCD spectrum that have been observed experimentally. 
A mode that is somewhat heavier, but not that much, than vector and axial vector modes is the 7Ti(1400) meson, with 
J PC — 1 h and a mass m wi ~ 1354 ± 25 MeV (l7j. A peculiarity of this mode is that it cannot be predicted within 
the valence quark model, or in other words a simple quark bilinear operator would not create this kind of mode. An 
operator with the right quantum numbers would involve also a gluon field ipFij^ip. Then, in order to include mesons 
with the quantum numbers of ttx, we would have to introduce a field dual to the dimension-five operators tpF^j^ip, 



action, since there is no self-duality constraint for these operators 



and ipFfj, u ip. The obvious candidate is again a complex two-form field, with bulk mass m l = 9 and no Chcrn-Simons 
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